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ABSTRACT 

We investigate the longevity of broad bridge features in position-velocity diagrams 
that appear as a result of cloud-cloud collisions. Broad bridges will have a finite 
lifetime due to the action of feedback, conversion of gas into stars and the timescale 
of the collision. We make a series of analytic arguments with which to estimate these 
lifetimes. Our simple analytic arguments suggest that for collisions between clouds 
larger than R ^ 10 pc the lifetime of the broad bridge is more likely to be determined 
by the lifetime of the collision rather than the radiative or wind feedback disruption 
timescale. However for smaller clouds feedback becomes much more effective. This is 
because the radiative feedback timescale scales with the ionizing flux Niy as liJ 
so a reduction in cloud size requires a relatively large decrease in ionising photons 
to maintain a given timescale. We find that our analytic arguments are consistent 
with new synthetic observations of numerical simulations of cloud-cloud collisions 
(including star formation and radiative feedback). We also argue that if the number 
of observable broad bridges remains ^ constant, then the disruption timescale must 
be roughly equivalent to the collision rate. If this is the case our analytic arguments 
also provide collision rate estimates, which we find are readily consistent with previous 
theoretical models at the scales they consider (clouds larger than about 10 pc) but are 
much higher for smaller clouds. 

Key words: stars: formation - ISM: kinematics and dynamics - ISM: clouds - ISM: 
Bubbles - galaxies: ISM - galaxies: star clusters 


1 INTRODUCTION 

Collisions between giant molecular clouds in galactic discs 
are a candidate mechanism for triggering the formation of 
massive stars and potentially even massive star clusters. 
When two clouds collide, a dense layer of material rapidly 
forms which might provide the conditions required for high 
mass star formation. This process has been studied in nu- 


merical simulations by, for example 

Stone (1970 

), Habe & 

Ohta (19921, 

Ricotti et al. (1997), Anathpindika ( 

2010), In- 

one h Fukui 

(2013), Takahira et al. (2014) and 

Wu et al. 

(2015), which show that collisions can indeed create the 


et al. (2015) used numerical simulations of gas evolving in a 


galactic potential to estimate a collision rate of about 1 per 
lOMyr for clouds larger than about 10 pc. Similar estimates 
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of the collision frequency for clouds larger than 10 pc have 
also been obtained by Tasker & Tan (2009), Tasker (2011) 


and Tasker et al. (2015). 


Observationally, many massive stars and star clusters 
have now been identihed as possibly resulting from cloud 
collisions. Two particularly compelling examples are M20 
( Torii et ar]|2011 ) and NGC3603 ( Fukui et al.||2014 ). These 
are star clusters which appear to be at a junction between 
two clouds with line of sight velocities differing by approx¬ 
imately 20km s“^. Of course it is possible that such a con- 
hguration could result from an independent star formation 
scenario and a chance alignment of two other clouds along 
the line of sight. However analysis of the (J=3-2, 2-1 
and 1-0) lines demonstrated that the components of these 
two separate clouds closest to the stellar cluster are being 
heated by the cluster - implying that they must be spatially 
coincident. 

Another signature of cloud-cloud collisions is a broad 
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Figure 1. A schematic of a collision-observer system and a car¬ 
toon p-v diagram showing a broad bridge feature. Different colour 
components in the collision schematic correspond to the different 
colours on the p-v diagram. 



Declination 


bridge feature in position-velocity diagrams (hereafter p-v 
diagrams). This is two velocity peaks, spatially coincident 
but separated in velocity, with intermediate intensity gas 
at intermediate velocities. In Haworth et al. (2015) we pro¬ 
duced synthetic p-v diagrams from a range of different sim¬ 
ulations of molecular clouds, including cloud collisions and 
isolated clouds with radiative feedback. We found that in¬ 
deed broad bridges appeared in our cloud collision models, 
but did not arise in any of the simulations of isolated clouds 
with radiative feedback (though in principle, a broad bridge 
could possibly arise given a favourable configuration of gas 
clouds merely coincident along the line of sight). [HawortlT 


et al. (2015) also compared the simulated p-v diagrams with 


observed broad bridge features towards M20, a site of colli¬ 
sion according to Torii et al. (2011). 

Broad bridges in p-v diagrams clearly can arise from 
cloud-cloud collisions, however it is currently unclear for 
how long they should survive. Feedback from winds and the 
ionising radiation field of any OB stars that form, as well as 
the conversion of gas into stars and the finite lifetime of the 
collision will presumably limit the amount of time that broad 
bridge features are visible for. In this paper we combine a 
series of analytic arguments with new synthetic observations 
of cloud collision models, both with and without radiative 
feedback, to try and address this question. 


2 ANALYTIC ARGUMENTS REGARDING 
THE TIMESCALE FOR DISRUPTION OF 
BROAD BRIDGES 

In this section we explore the time for which a broad bridge 
feature might survive. There are many processes that might 
act to remove the broad bridge such as the collision being 
completely braked or ending as one cloud punches through 
the other, radiative/mechanical feedback from massive stars 
and the conversion of gas into stars. 


Figure 2. An example CO (J=l —)■ 0) p-v diagram of a possible 
cloud-cloud collision towards M20 ( [Torii et al.||201l] ). This p- 
V diagram is produced using data taken using Mopra and was 
presented in Haworth et al (2015). 


2.1 Preamble: what is the broad bridge 

We begin by recapping what the broad bridge feature is in 
the context of a collision. Figure shows a schematic of a 
cloud-cloud collision, as well as a cartoon of the p-v dia¬ 
gram the observer would see. A smaller cloud (blue) drives 
into a larger cloud (red). The entirety of the smaller cloud 
undergoes collision quite quickly, resulting in a compressed 
layer that continues to move into the larger cloud. Between 
the bulk velocities of the compressed layer and larger cloud 
there is intermediate velocity gas (yellow). For the observer 
viewing angle in the schematic, the observer sees some red- 
shifted component moving away from them, a blue shifted 
component moving towards them and the intermediate ve¬ 
locity gas. In the p-v diagram this manifests itself as two 
peaks along the velocity dimension, separated by lower in¬ 
tensity intermediate velocity emission - the broad bridge 
feature. An example of a p-v diagram with a broad bridge 
towards in M20, taken with Mopra, is given in Figure 

The intermediate velocity gas is replenished as long as 
the collision is still occurring, so in order to remove the broad 
bridge feature the collision either needs to end (being fully 
braked or with one cloud punching right through the other), 
or the compressed dense layer resulting from the smaller 
cloud needs to be at least partially removed. 


2.2 Details of the collision 


We consider a collision between two clouds in a scenario sim¬ 
ilar to that discussed by Habe & Ohta ( |1992| . One cloud of 
radius Ri collides with a second, larger (R 2 ), cloud at veloc¬ 
ity Vc. The entirety of the small cloud has been compressed 
after a time approximately given by tc ~ Rijvc (this is a 
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lower limit, since it assumes that the post-collision veloc¬ 
ity is small). We work in a frame in which the larger cloud 
is stationary. After tc the compressed layer is still moving 
with some finite bulk velocity relative to the larger (in our 
frame, static) cloud, meaning that the broad bridge is still 
observable for some subsequent timescale (to be determined 
below). Our considerations here are only very approximate, 
we do not include the effects of turbulence or instabilities 


that can arise in a compressive flow (e.g. McLeod & Whit- 
worth|20T3 ). 


2.3 The collision timescales 

The broad bridge will disappear if the clouds brake to the 
extent that the are no longer separated in velocity by more 
than half the sum of their turbulent velocity dispersions. If 
this happens then only a single feature will appear in the 
p-v diagram. The broad bridge may also disappear if the 
smaller cloud punches right through the larger. To estimate 
these timescales we require a description of the smaller cloud 
velocity as a function of time. 


2.3.1 Momentum driven motion in the limit of eonstant 
deeeleration 


The smaller cloud (compressed layer) collides with the larger 
cloud in a frame in which the large cloud is initially static. 
We assume that all material from the larger cloud encoun¬ 
tered by the compressed layer is swept up to travel with it. 
By conservation of momentum, this sweeping up of material 
decreases the velocity in the compressed layer and will lead 
to the collision being completely stalled after a distance 


^ _ 2piRiVc 

P2v{t) ■ 


( 1 ) 


Assuming constant deceleration from the collision velocity 
Vc to zero 


dv _ Vc _ P2v{t) 

dR Dstag ^piRi 

which yields 

dv dv dR p 2 v{t)‘^ 

dt dR dt 2piRi 


( 2 ) 

(3) 


integrating this gives us a crude estimate of the velocity of 
the smaller cloud relative to the larger as a function of time 


v(t) = 


1 


P2t I J_ 

2piRi Vc 


2RipiVc 
p2tvc + 2piRi 


(4) 


This is shown for the 20 k ms ^ collision between large clouds 
(prep) as well as the 10kms“^ collision 


from Shima et al. 


model fromiTakahira et M (2014) in Figure]^ The Shima 


et al. ( prep| ) models include sink formation and radiative 
feedback whereas the |Takahira et al. ( 2014| models do not. 
Included in Figure are also points inferred from synthetic 
molecular line data cubes of the simulations (these synthetic 
observations will be discussed in more detail below). There 
are three time regions denoted on each plot. The leftmost re¬ 
gion (t < tc) is prior to the crossing timescale of the smaller 
cloud. In this region we don’t necessarily expect good agree¬ 
ment because only some small fraction of the smaller cloud 
is undergoing collision, which the analytic model does not 



Time, Myr 



Figure 3. The velocity of the compressed layer as a function 
of time estimated using momentum balance and by assuming 
constant deceleration. The points are extracted from p-v di¬ 
agrams/line profiles from synthetic data cubes of simulations 
(which we will discuss in section 3). The width of the secondary 
peak in the broad bridge, or ease with which two distinct veloc¬ 
ity components can be identified determines the size of the error 
bars. The left hand vertical line denotes the time after which the 
small cloud is compressed. The right hand vertical line denotes 
the time after which the smaller cloud has punched through the 
larger. 


account for. In the middle time region the whole small cloud 
is compressed into a layer travelling with a similar bulk ve¬ 
locity (ignoring turbulence), we therefore expect quite good 
agreement with our simple analytic estimate in this regime. 
The right hand time region is beyond the timescale that we 
predict the smaller cloud has exited the larger. In this re¬ 
gion it is again difficult to estimate the velocities from p-v 
maps since the broad bridge is completely removed or un¬ 
clear. The derived estimates in this final time zone are made 
using a combination of p-v maps, the raw simulation data 
and line profiles. 

The simple analytic approximation reproduces the val¬ 
ues seen in the simulations to within (at worst) ~20 per cent 
and will therefore be a suitable tool for giving a first order 
estimate of the timescale over which the broad bridge sur- 
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vives due to the overall collision lifetime. Although after tc 
the velocities predicted by our simple relation lie within the 
uncertainties (which are mostly based on how broad the two 
peaks of the broad bridge are) the simulations typically ex¬ 
hibit lower velocities than the analytic model predicts. This 
could be due to turbulent ram pressure braking the cloud 
further. The analytic solution also does not include gravity, 
which could affect the motion. 

Equation integrated until the collision velocity is less 
than half the sum of the two clouds turbulent velocity dis¬ 
persions 

v{t) < + Au 2 ), (5) 


or until the compressed layer escapes the larger cloud (tra¬ 
verses R 2 ) gives us our first broad bridge lifetime es¬ 

timate. For the latter criterion we obtain an exit timescale 
for the smaller cloud of 


te 


2piRi r .p2{Ri-\-R2) 

— - exp(^-^---^ 

p2Vc I 2piRi 



( 6 ) 


This is our first broad bridge lifetime estimate, we will utilise 
it in section |3] 


2.4 Radiative feedback disruption timescale 


During the collision material is compressed into a dense 
layer, which continues to travel into the larger cloud, giv¬ 
ing rise to the secondary feature in the p-v diagram. It is 
also within this dense region that star formation (possibly 
massive star formation) will occur. However, once stars have 
formed their ionising radiation will result in bubbles of hot 
gas that might disrupt the dense layer, potentially render¬ 
ing it unobservable (at least in molecular lines) once some 
fraction of it has been ionised. 

For a star emitting Niy ionising photons per second into 
a medium of number density n a bubble of ionised gas re¬ 
sults, with a “Stromgren” radius ( Stromgr'^|1939 ) 

^ / miy 


(7) 


where as is the hydrogen case B recombination coefficient 
(that for recombinations in to all states other than the 
ground). This bubble is hot (10^ K) and overpressured rel¬ 
ative to the ambient medium, meaning it expands. For our 
purposes, the extent of the H ll region (with sound speed Ci) 
under such so-called “D-type” expansion can be adequately 
described by the relation derived by |Spitz^ ( |1998 ) 

, 4/7 


r{t) = ' 


1 + 


7cit 

4rs 


( 8 ) 


(though the interested reader is directed towards Bisbas 
[eTallpoTsl l, for a discussion on the shortcomings and in¬ 
terpretation of the Spitzer solution, including how it differs 
from the solution derived by Hosokawa and Inutsuka, 2006). 
As the ionised region expands it will disrupt a larger volume 
of the material that makes up the broad bridge feature, until 
at some point it renders it undetectable. 

We assume that the compressed layer has a length ap¬ 
proximately the same as the smaller cloud diameter, 2 x Ri 
and that the width is less than the Stromgren radius (see 
appendix A). If this is the case, and we assume that the 


stars are centrally condensed in a cluster, then we only need 
to calculate the extent along the length of the layer that 
is ionised as a function of time. If the broad bridge is un¬ 
observable after some fraction / of the compressed layer is 
ionised, then we can calculate the broad bridge lifetime in 
the presence of radiative feedback using fRi = r(t), giving 




4rs 

7ci 



(9) 


We can recast in terms of collision variables. If the shock 
is isothermal, then the number density just behind the shock 
is 


n = 


P2 

IJlH 



( 10 ) 


where ul, P 2 and C 2 are the dense layer velocity and the 
larger cloud density and sound speed respectively. Substi¬ 
tuting equation 1 10| into equation results in 


yAnasp^vi J 


( 11 ) 


which combined with equation gives us a radiative feed¬ 
back disruption timescale as a function of the large cloud 
density and collision velocity 




( fRi ^ 

7/4 

4 

f SNiym^C2 \ ^ 

- 1 

7ci 

{AnanplvlJ 

_ \ \ 47 raBP 2 Rt ) ' 


( 12 ) 

At the first level of approximation, we might set the dense 
layer velocity equal to the collision velocity (vl = Vc). In¬ 
deed if these equations are applied to observations, we sug¬ 
gest just using the observed collision velocity. However, for 
our discussion here it is more accurate to use the mean veloc¬ 
ity over the lifetime of the broad bridge feature. To estimate 
the mean velocity, we first calculate an initial estimate of 
tr using the velocity at the onset of collision Vc- We then 
calculate the mean velocity over this initial tr timescale es¬ 
timate using equation Lastly, we iterate over tr and vl 
until the fractional change in mean velocity is less than 1 per 
cent. Typical mean velocities over the radiative timescale are 
around 40 to 60 per cent of the initial collision velocity. 

We plot the disruption timescale as a function of ion¬ 
ising flux for different cloud/collision parameters in Fig¬ 
ure 1^ each assuming that the clouds have a density of 
100 m/f cm“^. We assume that the ionised and neutral gases 
have temperatures (sound speeds) of 10^ and 10 K (11.5 
and 0.19 km/s) respectively. We also assume that the crit¬ 
ical fraction of the cloud length that needs to be ionised 
is / = 0.5. Different lines in Figure ^ correspond to differ¬ 
ent smaller cloud radii and collision velocities in the upper 
and lower panels respectively. Since the radiative disruption 
timescale scales as R^^ but only a reduction of the 

cloud size by a factor of 10 requires a reduction in the ion¬ 
izing luminosity by a factor of 10^ to maintain the same 
disruption timescale. We therefore expect broad bridges to 
be disrupted more readily by radiative feedback in collision 
between smaller clouds, even if the ionizing luminosity is 
substantially lower. We will compare these relations with 
simulations in section]^ 
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Figure 4. The timescale for a broad bridge to be disrupted by ra¬ 
diative feedback as a function of ionising flux. The upper panel is 
for a 20kms“^ collision for different cloud radii. The lower panel 
is for a 10 pc cloud and different collision velocities. These plots 
assume that both colliding clouds have a density of 100 mn cm“^. 



Figure 5. The timescale for a broad bridge to be disrupted by 
stellar wind feedback as a function of wind luminosity for different 
cloud radii (at 20km s“^ collision, top panel) and collision veloc¬ 
ities (for a cloud of 10pc radius, lower panel). The plots assume 
that both clouds have a density of 100 mn cm“^. 


2.5 Stellar wind feedback disruption timescale 


We can estimate the effect of stellar wind feedback using 
a similar approach to that for radiative feedback. That is 
we can compare the size of a wind blown bubble with the 
size of the compressed layer. We base our treatment of the 
expanding wind-blown region on the expansion in a cold 


uniform medium presented by Ngoumou et al. (2015), which 


in turn is based on Steigman et al. (1975) and Lamers & 


Cassinelli (1999). In this model the bubble radius is given 


by 




(13) 


where Vw and ps are the wind mass loss rate, velocity 
and slab density (all in cgs units) respectively. Comparing 
equation with some critical fraction of the layer size fRi, 
and assuming an isothermal shock, gives a time at which 
disruption of the broad bridge due to winds takes place of 


tw - 


fRlpl^^VL 

0.6889M„^^%i/^C2 


(14) 


Note that this assumes that the total wind luminosity 
goes in to disrupting the broad bridge. 


Plots of the broad bridge disruption timescale as a func¬ 
tion of wind luminosity for a variety of cloud sizes and colli¬ 
sion velocities are given in Figure Note that in these plots 
we assume a wind velocity of 2000 km s“^ and a gas density 
of lOOmHcm”^. We again calculate the mean dense layer 
velocity using the iterative procedure discussed in section 

EH 

By comparing equation with equation for systems 
where [fRi/rsY^^ ^ 1 our radiative feedback timescale 
dominates over that for winds when 




10 ~^Mq yr“ 


< 1.7 



Nly 

1048S-1 


1/2 


(15) 


where our assumed values are that the wind velocity is 
2000 km s“^, the ionised and neutral gas temperatures are 
10^ and 10 K respectively, / = 0.5 and case B recombina¬ 
tion coefficient of 2.7 x 10“^^ cm^ s“^. According to previous 


numerical simulations by Dale et al. (2013) and Dale et al. 


(2014) radiative feedback is typically expected to dominate 
over winds. 
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2.6 The star formation disruption timescale 

The gas mass depletion through star formation is roughly 
given by the time integral of the average star formation rate 
M. The broad bridge will be disrupted when some fraction F 
(which is probably different from the fraction of the length 
/ discussed in the feedback timescales) of the dense layer 
mass Ml is converted into stars 

jMdt = FMi. (16) 

then the timescale for star formation to remove the broad 
bridge is 

t,f = FMi/M. (17) 

Unfortunately to, F, Mi and M are all rather uncertain, 
so a more useful relation is difficult to derive. In subsequent 
comparison with numerical models in this paper we therefore 
extract the values of to and M from the simulations. We 
assume that Mi is the smaller cloud mass and T is 0.5. This 
is the most uncertain timescale in this paper. 


2.7 Weaknesses of the analytic arguments 


Our arguments made above are purely intended to provide 
rough lifetime estimates and to build an understanding of 
how feedback might disrupt the broad bridge feature. They 
may also assist in a basic interpretation of real observations. 
However, they have a number of shortcomings 


(i) Simulations by Dale & Bonnell (2011) find that ra¬ 
diative feedback can be ineffective if the H ll region is sub¬ 
stantially matter bounded. That is, if photons are able to 
stream freely from the ionised bubble. This would make ra¬ 
diative feedback much less effective than we predict and so 
will increase tr- 

(ii) Regarding both winds and ionising radiation, we as¬ 
sume that the sources responsible for feedback are co¬ 
located. In reality they will be distributed throughout the 
colliding layer. 

(hi) At present there is no consideration of supernovae, 
which might also be effective at removing broad bridges after 
about 5Myr. 

(iv) We do not account for the effects of turbulence and 
assume constant deceleration in our collision timescale esti¬ 
mate. 

(v) We assume that the collision is isothermal. In real¬ 
ity the compressed layer density will likely larger than that 
given by the isothermal shock condition (equation 10), at 
least in the absence of turbulence. 


3 NUMERICAL COMPARISON 

We now study the evolution of the broad bridge feature in 
two numerical models and compare with our expectations 
from the analytic expressions. One of the models is without 
any form of feedback or sink formation, the second includes 
radiative feedback and sink formation. Neither include feed¬ 
back from winds. 


Table 1. Summary of the parameters of the feedback model from 
Shima et al (in prep). 


Parameter 

Value 

Description 

Ri 

14.45pc 

Smaller cloud radius 

R2 

27.93pc 

Larger cloud radius 

Pi 

24.4 niH cm“^ 

Smaller cloud density 

P2 

11.79 ruH cm“^ 

Larger cloud density 

Vc 

20 km 

Collision velocity 

Avi 

1 km 

Smaller cloud turbulent velocity 

Av2 

1.49 km 

Larger cloud turbulent velocity 

Q 

3 X 

Ionizing flux per solar mass 

Qtot 

3 X 

Total ionizing flux after 2.4 Myr 
of collision 


Table 2. Summary of the parameters of the feedback model from 
Takahira et al (2014). 


Parameter 

Value 

Description 

Ri 

3.5pc 

Smaller cloud radius 

R2 

7.2pc 

Larger cloud radius 

pi 

47.4 mi(f cm“^ 

Smaller cloud density 

P2 

25.3 m/f cm“^ 

Larger cloud density 

Vc 

lOkms-i 

Collision velocity 

Avi 

1.25 kms“^ 

Smaller cloud turbulent velocity 

Av2 

1.71 kms“^ 

larger cloud turbulent velocity 


3.1 Numerical method 

3.1.1 Hydrodynamic models 


The hydrodynamic models that we post process are simula¬ 
tions of cloud-cloud collisions by Takahira et al. (2014) and 
Shima et al (in prep). These models were both run using 


the grid based radiation hydrodynamics code ENZO (Bryan 
et al.|2014| |. They both consider head-on collisions between 


two non-identical, turbulent clouds. 

The collision with feedback and sink particles is the 
20kms“^ collision model from Shima et al (in prep). A 
summary of the cloud parameters is given in Table Once 
sink particles form in this simulation they are assumed to 
have an ionizing luminosity of 3 x lO^^photonss”^ 
this model both clouds are travelling in opposite directions 
at 10kms“^. We consider snapshots at 0.4, 2, 2.4, 4.4 and 
5.4 Myr after the onset of collision. At 2.4 Myr there is ap¬ 
proximately 100 Mq in sink particles, implying an ionizing 
flux of 3 X lO^^photonss this is the value we use in our 
radiative feedback disruption timescale estimate. 

The collision without feedback or sink particles (i.e. 
without star formation) is the 10 km s“^ collision model from 
Takahira et al.| ( |2014| . A summary of the clouds parameters 


IS given m Table In this model the collision is in such a 
frame that the larger cloud is stationary. We consider snap¬ 
shots at 0.1, 0.2, 1, 1.6, 2.6 and 4.4 Myr after the onset of 
the collision. 


3.1.2 Radiative transfer 

As for our original synthetic observations in [Haworth et ah] 
(2015) we use the TORUS radiation transport and hydrody¬ 


namics code to post-process simulations in this paper (Har- 
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Figure 6. Various timescales for associated with the broad bridge 
for the 20kms“^ collision model from Shima et al (in prep), tgf, 
tr and tw are the time it takes for star formation, radiative feed¬ 
back and wind feedback to quench the broad bridge, te is the time 
it takes for the small cloud to punch through the larger cloud. 


Wind mass loss rate, Mg^/yr 

10 '^ 10 '® 10 '® 10 "'' 10 ‘® 



Figure 8. Various timescales for associated with the broad bridge 
for the 10kms“^ collision model from Takahira et al (2014). tgf, 
tr and tw are the time it takes for star formation, radiative feed¬ 
back and wind feedback to quench the broad bridge, te is the time 
it takes for the small cloud to punch through the larger cloud. 


ries|20QQ] Kurosawa et al.|2004| Haworth Harries|2012 ). In 


particular we use the molecular line transfer algorithm dis¬ 


cussed by Rundle et al. (2010), Haworth et al. (2013). TORUS 


is capable of non-LTE statistical equilibrium calculations, 
which were used in Haworth et al. (2015), however since we 


are prostprocessing a relatively large number of snapshots 
in order to study the time evolution of broad bridges we 
assume LTE when calculating the molecular level popula¬ 
tions in this paper. This does not alter the qualitative p-v 
diagram (we are only studying the morphological evolution) 
and substantially reduces our computational expense. With 
the level populations calculated, a synthetic data cube is 
generated using ray tracing ( Rundle et al.[|2010 ). This data 
cube is manipulated using the STARLINK image analysis tool 
GAIA to produce p-v maps. In this paper we produce syn¬ 
thetic (J=l ^ 0) data assuming an abundance relative 

to hydrogen of 8 x 10“^. 


3.2 Results of synthetic observations 

3.2.1 20km s~^ collision model from Shima et al (in prep) 


We first apply our analytic arguments to a 20kms“^ col¬ 
lision between two large clouds studied by Shima et al (in 
prep). The parameters of this model are given in Table 
In Eigurej^we plot the broad bridge disruption timescales 
based on the exit timescale (equation |^, wind/radiative 
feedback (equations and 14) and and star formation rate 
(equation pT| for this model. We allow the star forma¬ 
tion rate and wind mass loss rate to vary and use the 
known ionizing flux of stars in the simulation at 2.4 Myr 
(3x10^® photons s“^). Erom these disruption timescales we 
expect the broad bridge to be first disrupted after about 
3.5 Myr because the smaller cloud has exited the larger (see 
the red small-dotted horizontal line in Eigure[^. High star 
formation rates or wind mass loss rates are required to dis¬ 
rupt the broad bridge on a comparable timescale. We expect 


that radiative feedback would require about 11 Myr to clear 
the broad bridge. 

Eigurej^shows the p-v diagram for this model at 0.4, 2, 
2.4 and 4.4 Myr after the onset of collision. In the upper left 
panel, the clouds have only just begun colliding. This snap¬ 
shot is at a time earlier than the timescale for crossing the 
small cloud at the collision velocity. At this stage in the col¬ 
lision the p-v diagram just looks like clouds coincident along 
the line of sight (though there is some very weak intermedi¬ 
ate velocity emission). As the collision progresses the broad 
bridge forms, getting harder to identify as the two intensity 
peaks approach one another and as the smaller cloud’s sig¬ 
nature becomes more disordered. By 4.4 Myr (bottom right 
panel) the broad bridge is obviously not present, consistent 
with the result expected from our simple analytic estimates. 
Although the broad bridge is not observable, there is dis¬ 
ordered emission from gas over a 10kms“^ velocity range 
at this time. By 4.4 Myr the larger cloud’s feature has also 
started to be significantly disrupted. The larger feature also 
appears to decrease in velocity dispersion with time (it’s 
thickness along the velocity axis decreases). This is due to 
dissipation of the turbulence, either by collisions or dilution 
in the wider ISM. 

Our analytic arguments give a broad bridge lifetime 
consistent with the observations to within at least 1 Myr 
in this case. According to these arguments the broad bridge 
lifetime in this scenario is determined by the collision life¬ 
time, rather than feedback. Though it is likely that conver¬ 
sion of gas into stars, feedback and the collision lifetime do 
all cumulatively contribute, meaning that the collision life¬ 
time estimate is likely an upper limit on the broad bridge 
lifetime. 


3.2.2 10km s ^ collision model from Takahira et al (2014) 

We now apply our analytic arguments to the 10kms“^ col¬ 
lision model from Takahira et al. (2014), the parameters 
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Figure 7. (J=l —)■ 0) synthetic p-v diagrams of snapshots from the 20kms ^ collision model from Shima et al (in prep). The 

panels are for snapshots at 0.4 (top left), 2 (top right), 2.4 (bottom left) and 4.4Myr (bottom right). 


of which are summarised in Table Unlike in the previ¬ 
ous model, this does not include any radiative feedback or 
star formation. The disruption timescales for this model are 
shown in Figure where we assume an ionizing luminosity 
an order of magnitude lower than that in the previous model 
(10^^ photons/s). In this collision between smaller clouds our 
analytic arguments predict that star formation as well as 
radiative and wind feedback will be much more effective at 
disrupting the broad bridge - to the extent that they can 
readily dominate over the collision timescale. Recall that 
since the radiative disruption timescale scales as but 

only ^ a reduction of the cloud size by a factor of 10 

requires a reduction in the ionizing luminosity by a factor of 
10^ to maintain the same disruption timescale. The simula¬ 
tions in Takahira et al. (2014) did not include star formation 


or feedback, so we expect the broad bridge to survive un¬ 
til the escape time. Since the predicted escape timescale is 
smaller than the radiative one we do not expect radiative 
feedback would have been important to the broad bridge 
lifetime had it been included. 


P-v diagrams for snapshots of this model at 0.1, 1, 1.6 


and 2.6 Myr are shown in Figure There is similar evolu¬ 
tion to the other collisional model, only with less small scale 
clumpy structure due to the lack of conversion of gas into 
stars, as well as a lack of feedback. The bottom right panel 
of Figure is 0.2 Myr after the exit timescale predicts the 
broad bridge should be disrupted. Although some remnant 
of a secondary intensity peak remains it has been predom¬ 
inantly disrupted and is no longer so “broad”. The actual 
exit time in the simulation is about 2.8 Myr, compared to 
our analytic estimate of 2.2 Myr. This difference is consistent 
with our findings in section 2.3.1, where the actual collision 
velocity in the simulations was typically slightly lower than 
the analytic relations predict. This is due to lack of treat¬ 
ment of processes such as turbulent pressure and gravity. 


4 DISCUSSION 

4.1 On the cloud—cloud collision rate 

We have developed tools for estimating broad bridge life¬ 
times which, where testable, seem to be consistent with 
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Position, pc 


Position, pc 


Figure 9. (J=l —)► 0) synthetic p-v diagrams of snapshots from the lOkms ^ collision model from Takahira et al (2014). The 

panels are for snapshots at 0.1 (top left), 1 (top right), 1.6 (bottom left) and 2.6Myr (bottom right).. 


our numerical models. In general we expect broad bridge 
lifetimes of around 1-5 Myr. To date, seven broad bridge 
features have been identified in the Galaxy. Six of these 


are strong: M20 (Torii et al. 2011 Haworth et al.| 2015), 
RCW120 UTorii et al.||2015|l, NGC3603 (|Fukui et al.pOWl, 
RCW38 (|E\ikui et al.|2015||, NGC6334 and NGC 6357. West- 


erlund 2 also shows a weak broad bridge feature (Furukawa 
et al.||2009 ). Dobbs et al. (2015) predict collisions on of or¬ 
der one every 10 Myr for clouds of size 10 pc or larger, which 
would be inconsistent with the amount of broad bridge fea¬ 
tures currently observed (which are a lower limit, given that 
more may yet be observed) if our typical broad bridge life¬ 
times are correct. 

However, in order for the number of broad bridges to not 
be extremely large (or zero), there should always be approx¬ 
imately an equilibrium number of broad bridges. In order 
to achieve a steady state, the rate at which broad bridges 
are formed must be approximately equal to the destruction 
rate. We can hence use our analytic arguments to estimate 
the equilibrium collision rate as a function of cloud size. 
We do this in Figure 10 for 10kms“^ collisions. We assume 
that the larger cloud radius is twice that of the smaller. We 


assume that the smaller cloud has a density that scales as 

n = max(-2.1* T_|_ 54.75^ 10.) (18) 

pc 

(which is based on the radii and densities of the clouds in 
the dynamical models discussed here) and that the larger 
cloud is half of this density. We also assume that the ionizing 
luminosity of the resulting stars varies with colliding cloud 




3 X 

14.45 pc 


(19) 


which is motivated by the Shima et al (in prep) simulation. 
We only consider the maximum required collision rate based 
on the exit timescale (equation]^ and the radiative feedback 
timescale (equation . 

For clouds larger than about 5 pc the exit timescales 
require collision rates consistent with those calculated by 


Dobbs et al. (2015), but for smaller clouds (not identified in 
the Dobbs et al.||2Ql^ estimates) the collision rates are ex¬ 
pected to be much higher (this is required by both the exit 
and radiative feedback timescales). Of the systems observed 
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Figure 10. The time between cloud-cloud collisions as a function 
of cloud radius assuming that an equilibrium number of broad 
bridge features is maintained. This is for collision parameters dis¬ 
cussed in section 4.1. 


to host broad bridges, NGC6334 and NGC6357 are approx¬ 
imately 10 pc in size, however all others are less than 4 pc. 
For the example of M20 which is about 2 pc in size, from 
our equilibrium consideration we expect about 4 such colli¬ 
sions per Myr. This is in contrast to the Dobbs et al. (2015) 
general value of about 1 per 10 Myr. It is therefore possible 
that high cloud-cloud collision rates, leading to massive star 
formation, are present in the Galaxy that are not contradic¬ 
tory to those rates expected from Dobbs et al. (2015). Higher 
collision rates for smaller clouds could be tested using new 
galactic models that are able to identify clouds below 10 pc 
in size. If higher collision rates for smaller clouds are found, 
this would call for further detailed simulations of the colli- 
sional process. In particular, to understand the link between 
the properties of the colliding clouds (e.g. radius, density, 
turbulence, collision velocity) and any subsequent star for¬ 
mation (e.g. star formation efficiency and mass spectrum), 
which might be compared with observations of candidate 
collision sites. 


5 SUMMARY AND CONCLUSIONS 

We study the lifetimes of broad bridge features in p-v 
diagrams. We do so using a combination of analytic argu¬ 
ments and by post processing hydrodynamic and radiation 
hydrodynamic simulations of cloud-cloud collisions. We 
draw the following main conclusions from this work. 

1) We derived simple analytical tools to estimate broad 
bridge lifetimes. The lifetime of the broad bridge is con¬ 
strained by either the collision timescale, the timescales for 
disruption by winds or radiative feedback and the timescale 
for the depletion of gas through star formation. 

2) The radiative feedback lifetime tn scales with the 

ionising flux as Ir oc and the smaller cloud radius 

as tn oc rV^. This implies that to maintain a given 


radiative feedback disruption timescale for smaller clouds, 
the ionizing flux needs to be more drastically reduced. Due 
to this effect, although larger clouds have their broad bridge 
lifetimes determined by the collision timescale, for smaller 
clouds feedback can play a dominant role in determining 
the broad bridge lifetime. 


3) If we assume that there should be some constant number 
of broad bridges with time, then the disruption timescales 
should be roughly equal to the time between collisions at a 
given cloud size. If this is the case our disruption timescales 
give collision rate estimates. We find that such estimates 
are readily consistent with the 1 per 10 Myr collision rate 
estimates by [Dobbs et al. (2015) for clouds larger than 


about 6 pc. However for smaller clouds (which Dobbs et al. 


2015 


do not consider) we find much high collision rates, for 


example 4 per Myr for 2 pc radius clouds. 


4) Given conclusion 3 and the fact that most high mass 
star formation sites where broad bridges have been iden¬ 
tified are less than 4 pc in size, we suggest that a large 
number of smaller scale cloud-cloud collisions could be 
triggering a substantial amount of high mass star formation. 
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We can approximate the dense layer as a self-gravitating 
isothermal slab, for which the density distribution is given 
by 


Pl(z) = pL{0)sedd 
where the scale height is 
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(see e.g. Umekawa et ah 1999| Glarke fc GarsweTI||2007| ). 
We assume that the slab sound speed Cs is the same as 
for all neutral gas (i.e. the same as that in the neutral col¬ 
liding clouds). Equating this scale height to the Stromgren 
radius (equation [D we find that the assumption that the 
Stromgren radius is larger than the slab width applies wher- 
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Which, inserting standard values, can be expressed as 
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Where we remind the reader that is the dense layer veloc¬ 
ity and 77-2 is the density in the larger cloud. Eor the example 
of the Shima et al (in prep) collision (discussed here in sec¬ 
tion 3.2.1) we find a predicted slab scale height of 0.05 pc 
(so a slab width of about 0.1 pc) and that equation A4 is 


satisfied. We remind the reader that this analysis assumes 
that the ionising sources are centrally located. 

A similar consideration can be made for winds by com¬ 
paring equations |A2| and |13| which tells us that the wind 
will break out of the slab after a time 
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This paper has been typeset from a TgX/ DTeX file prepared 
by the author. 


APPENDIX A: THE THICKNESS OF THE 
DENSE LAYER 

Here we discuss the assumption that the slab thickness is 
less than the Stromgren radius. We begin by noting that 
this criterion is not essential and that so long as the layer is 
not much larger than the Stromgren radius (meaning that 
the Hii region will expand beyond the slab thickness quite 
quickly) our results are unaffected. 




















